Multiplicity scaling in ideal and viscous hydrodynamics 
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Using numerical results from ideal and viscous relativistic hydrodynamic simulations with three 
different equations of state, for Au+Au and Cu+Cu collisions at different centralities and initial 
energy densities, we explore the dependence of the eccentricity-scaled elliptic flow, V2/e, and the 
produced entropy fraction, AS /So, on the final charged hadron multiplicity density dNch/dy per 
unit transverse overlap area 5, {1/ S)dNch/dy. The viscous hydrodynamic simulations are performed 
with two different versions of the Israel-Stewart kinetic evolution equations, and in each case we 
investigate the dependence of the physical observables on the kinetic relaxation time. We find 
approximate scaling of V2/e and AS /So with {l/S)dNch/dy, with scaling functions that depend on 
the EOS and, in particular, on the value of the specific shear viscosity rj/s. Small scaling violations 
are seen even in ideal hydrodynamics, caused by a breaking of the scale invariance of ideal fluid 
dynamics by the freeze-out condition. Viscous hydrodynamics shows somewhat larger scale-breaking 
effects that increase with increasing tj/s and decreasing system size and initial energy density. We 
propose to use precision studies of these scaling violations to help constrain the shear viscosity rj/s 
of the quark-gluon plasma created in relativistic heavy ion collisions. 

PACS numbers: 25.75.-q, 12.38. Mh, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Recent numerical studies of relativistic hydrodynam- 
ics for dissipative fluids [1~4] have confirmed earlier es- 
timates [5, 6] that the "elliptic" anisotropic collective 
flow observed in non-central heavy-ion collisions is very 
sensitive to the shear viscosity of the matter formed in 
such collisions. Since ideal fluid dynamics (i.e. the as- 
sumption that viscosity can be neglected) provides a phc- 
nomenologically quite successful description of much of 
the soft hadron data collected from Au-|-Au collisions at 
the Relativistic Heavy Ion Collider (RHIC) [7-9], this 
implies strong constraints on the shear viscosity to en- 
tropy ratio rj/s [6, 10] and the thermalization time scale 
[11, 12] of the matter in the collision fireball. The con- 
clusion is that the quark-gluon plasma (QGP) created 
at RHIC is a strongly coupled plasma with almost per- 
fect liquid collective behavior [13-15] whose specific shear 
viscosity is lower than that of any previously known real 
fluid and consistent with a postulated lower bound of 
77/s > h/^ATrks) derived from the study of infinitely 
strongly coupled conformal field theories [16, 17] using 
the AdS/CFT correspondence and corroborated by ear- 
lier quantum mechanical arguments based on the uncer- 
tainty relation [18]. 

On the other hand, heavy-ion data at RHIC, SPS 
and AGS also show that ideal hydrodynamics gradually 
breaks down at larger impact parameters, for smaller col- 
lision systems, at lower collision energies, and away from 
midrapidity (see [9] for a review) . Much of this can be at- 
tributed to strong viscous effects in the late collision stage 



after the QGP has converted to hadrons [19]. However, 
due to uncertainties in the initial conditions for the fire- 
ball deformation, there is some room left for a non-zero 
value of the QGP viscosity [19]. To study this further 
requires a viscous hydrodynamic approach because the 
tool used to describe the non-equilibrium late hadronic 
dynamics (a classical cascade) is not well suited for the 
rapidly evolving, very dense matter in the early collision 
stages 

The motivation for the present paper is provided by 
the well-known systematic comparison of Voloshin et al. 
[20, 21] of elliptic flow data with ideal fluid dynamical 
predictions which suggests that the elliptic flow param- 
eter V2 scaled by the initial source eccentricity e, V2/S, 
while strongly deviating from ideal hydrodynamics at low 
multiplicities, still scales with the final multiplicity per 
unit overlap area: 



V2_ 1 dNch 
e S dy 



(1) 
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For ideal fiuids the right hand side is a direct measure 
of the initial entropy density [22]. The scaling (1) im- 
plies that all dependence on impact parameter, collision 
energy and system size can be, to good approximation, 
absorbed by simply taking into account how these control 
parameters change the final hadron multiplicity density. 
We will call this observation simply "multiplicity scaling 
of the elliptic flow" , where "elliptic flow" is a shorthand 
for the eccentricity-scaled elliptic flow V2/S and "multi- 
plicity" stands for {l/S)dNch/dy. 

Such a scaling is expected for ideal fluid dynamics 
whose equations of motion are scale invariant and where 
the eccentricity-scaled elliptic flow is therefore predicted 
[23, 24] to depend only on the squared speed of sound, 
= which describes the stiffness of the equation 



2 



of state (EOS) or "pushing power" of the hydrodynami- 
cally expanding matter. It has been known, however, for 
many years [25] that this ideal-fluid scaling is broken by 
the final freeze-out of the matter: if hadron freeze-out is 
controlled by hadronic cross sections (mean free paths) 
or simply parametrized by a critical decoupling energy 
density Cdec or temperature Idee, this introduces and ad- 
ditional scale into the problem which is independent of 
(or at least not directly related to) the initial geometry 
of the fireball. This breaks the above argument based 
on scale invariance of the ideal fluid equations of motion. 
We will show here that this also leads to a breaking of 
the multiplicity scaling of vije not only in the most pe- 
ripheral or lowest energy collisions, where freeze-out ob- 
viously cuts the hydrodynamic evolution short since the 
freeze-out density is reached before the flow anisotropy 
can fully build up [25], but even in the most central col- 
lisions at RHIC where freeze-out still terminates the hy- 
drodynamic evolution before the elliptic flow can fully 
saturate (see also [19]). 

The more interesting aspect of the experimentally ob- 
served scaling is, however, its apparent validity in regions 
where ideal fluid dynamics does not work (these encom- 
pass most of the available data [20]). Many years ago, 
simple scaling laws for the centrality dependence of ellip- 
tic flow were derived from kinetic theory in the dilute gas 
limit, where the particles in the medium suffer at most 
one rescattering before decoupling [26, 27]; these can be 
reinterpreted in terms of multiplicity scaling for vije. 
The dilute gas limit is expected to hold for very small 
collision systems, very large impact parameters or very 
low collision energies. More recently, a successful attempt 
was made to phenomenologically connect the dilute gas 
and hydrodynamic limits with a 1-parameter fit involv- 
ing the Knudsen number [28] . This fit works very well for 
Au+Au and Cu+Cu data from RHIC, but predicts that 
even in the most central Au+Au collisions at RHIC the 
ideal fluid dynamical limit has not yet been reached and 
is missed by at least 25% [28] . In the present paper we use 
viscous relativistic hydrodynamics to explore the multi- 
plicity scaling of vije in the phenomenologically relevant 
region. We conclude (not surprisingly since much of the 
available data is from regions where the viscous hadronic 
phase plays a large role [19]) that the multiplicity scaling 
data [20, 21] require significant shear viscosity for the 
medium, but also that viscous hydrodynamics predicts 
subtle scaling violations which seem to be qualitatively 
consistent with trends seen in the data (even if the ex- 
perimental evidence for scaling violations is presently not 
statistically robust) and whose magnitude is sensitive to 
the specific shear viscosity 77/s. This gives hope that fu- 
ture more precise data can help constrain the QGP shear 
viscosity through exactly such scaling violations. 

We should caution the reader that, similar to Ref. [28] 
which used a constant (time-independent) Knudsen num- 
ber, our viscous hydrodynamic calculations are done with 
a constant (temperature-independent) specific shear vis- 
cosity 77/ s. Neither assumption is realistic, and we expect 



77/ s in particular to show strong temperature dependence 
near (the critical temperature for the quark-hadron 
phase transition) and emerge from the phase transition 
with much larger values than in the QGP phase. Com- 
parisons between the results presented here and experi- 
mental data are therefore, at best, indicative of qualita- 
tive trends, and improved calculations, which in partic- 
ular match viscous hydrodynamics to a realistic hadron 
cascade below Tc, are required before an extraction of 
77/s from experimental data can be attempted. 

II. DISSIPATIVE FLUID DYNAMICS 

In this section we briefly review the viscous hydrody- 
namic equations that we solve, focussing on some differ- 
ences in the formulations used in previously published 
papers [1-4] which we investigate here further. For tech- 
nical details we ask the reader to consult these earlier 
papers. 

We focus on systems with exact longitudinal boost- 
invariance and use the code VISH2+1 [3] to solve numeri- 
cally for the expansion in the two dimensions transverse 
to the beam direction. As in [1-4] we consider only shear 
viscosity, neglecting bulk viscosity and heat conduction. 
(Bulk viscosity may become large near the QCD phase 
transition [29-31] and should thus be included in the fu- 
ture before comparing viscous hydrodynamics with ex- 
perimental data.) VISH2+1 solves the conservation laws 
for energy and momentum, dmT™" = (where dm is the 
covariant derivative in our curvilinear {T^x,y,r]) coordi- 
nate system [32]), with the decomposition 

T™" = eu"u" - pA"" + tt"". A'"" = (2) 

together with evolution equations for the viscous shear 
pressure tensor components tt"™: 

UTT = (TT —ZTja ) — (U IT +U IT jJJUk 

(3) 

Here D = u™'dm is the time derivative in the local co- 
moving frame, V™ = A'"'d; is the spatial gradient in that 
frame, and cr™" = V<™ii"> = i(V™w"+V"w™) - ^A""'e 
(with the scalar expansion rate = dkU^ = VfcU*'') is the 
symmetric and traceless velocity shear tensor, ujmn = 
V,iM„i — VmUn is the vorticity tensor, and a'-^b^^ = 
^{a"''b" + a"b"^) denotes symmetrization. Even though 
several components of the symmetric shear pressure ten- 
sor TT™" are redundant [32] on account of its traceless- 
ness and transversality to the flow velocity u™, VISH2+1 
propagates all 7 non-zero components and uses the trace- 
lessness and transversality conditions as checks of the nu- 
merical accuracy [3]. 

For a conformally symmetric fluid such as a massless 
quark-gluon gas, the temperature T is the only scale in 
the problem and therefore 77 ~ s T'^ and ~ ^/T, 
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hence j^T/t^ T^. In this Umit the first term in the 
second line of Eq. (3) can be written as [33] 

"r""^'^'' (^"') = +^^""(5^(1^^^) - 0) ■ (4) 

Tliis is tlic form used in Ref. [1]. 

It has recently been argued [34-37] that the r.h.s. of 
Eq. (3) should contain even more terms, at least for con- 
formal fluids in the strong coupling limit. Wc will not 
pursue this possibility here. 

Equations (3) are known as "Israel-Stewart (I-S) equa- 
tions" and based on an expansion of the entropy pro- 
duction rate to second order [33, 38-41] (macroscopic 
approach) or, in a microscopic kinetic approach using 
Grad's 14-moment method, of the phase-space distribu- 
tion function to first order in small deviations from local 
thermal equilibrium [38, 41, 42] (see also [43]). By intro- 
ducing a finite and sufficiently large relaxation time r^r 
for the evolution of the shear pressure tensor towards its 
Navier-Stokcs limit tt™" = 2rja™\ these equations elim- 
inate problems with acausal signal propagation at short 
wavelengths and the resulting numerical instabilities that 
famously plague the relativistic Navier-Stokes equation. 
A somewhat different approach to solving these problems 
was developed by Ottinger and Grmcla (0-G) [44] and 
has been used in [4] ; because a comparison of results ob- 
tained with the I-S and 0-G equations is non-trivial, we 
will leave that for a later study. 

Refs. [1-3] use different versions of Eqs. (3). P. & U. 
Romatschke [1] use the full set of terms displayed in (3) 
which we label as "full I-S equation". The last term in 
the second line of Eq. (3) involving the vorticity cannot 
be obtained from the macroscopic approach [32, 40] since 
it does not contribute to entropy production, but it fol- 
lows from the microscopic kinetic approach [38, 42]. For 
longitudinally boost-invariant systems with initially van- 
ishing transverse flow it is zero initially and was found in 
[1] to remain tiny throughout the fireball evolution. We 
can therefore remove it from consideration when compar- 
ing published results from the different approaches. The 
first term in the second line of (3) arises in this form from 
the macroscopic approach (2""^ order entropy production 
[32, 40]) but was neglected in our previous work [3], fol- 
lowing an argument in [32] that it is of second order in 
small quantities and therefore subdominant compared to 
the first two terms on the r.h.s. of Eq. (3). A simi- 
lar argument can be made for the last term in line 1 of 
(3) (which does not contribute to entropy production ei- 
ther), but Baier et al. [42] pointed out that this term is 
needed to preserve the transversality of tt™" during ki- 
netic evolution. In Ref. [3] we therefore kept all terms in 
the first line of Eq. (3) but dropped those in the second 
line; we call this here the "simplified I-S equation" . In 
Ref. [34] Baier et al. argued that for a conformally in- 
variant medium, such as a classical massless quark-gluon 
gas, the first term in the second line is needed to preserve 
the conformal invariance of the kinetic evolution equation 
and hence should not be dropped. We now understand 



that the arguments presented in [32] to neglect all but the 
first term on the right hand side of Eq. (3) were at best 
superficial since this term involves the difference between 
two first-order quantities and thus presumably needs to 
be counted as small of second order. 

Chaudhuri followed in his work [2] the approach ad- 
vertised in [32]; as a result, in his procedure tt™" must 
be expected to evolve away from transversality. He cir- 
cumvents this problem by evolving only the three lin- 
early independent components of tt™" and computing 
the rest from the tracelessness and transversality con- 
ditions [2]. The problem resurfaces, however, since now 
the results for all components of tt™" must be expected to 
depend on the choice of independent components which 
are evolved dynamically with the truncated equation (3). 
While these questions await quantitative study we note 
that Chaudhuri's results [2] appear to differ significantly 
from ours [3] . 

In the present paper we show many comparisons be- 
tween solutions obtained by using the "full I-S equations" 
with those from the "simplified I-S equations" . Although 
at sufficiently long wavelengths both have to agree in the 
Navier-Stokes limit r,r — * (up to issues of numerical 
stability), as inspection of Eq. (3) readily shows, they 
differ for non-zero r^r and will be seen to exhibit differ- 
ent degrees of sensitivity to r^r. This is of phenomeno- 
logical importance since r^r for the QGP is not known, 
and a strong sensitivity to this unknown parameter will 
negatively impact our ability to extract the QGP shear 
viscosity 77/s from experimental data. 



III. INITIAL CONDITIONS, FREEZE-OUT, 
AND EQUATION OF STATE 

For the present study, we initialize the expanding fire- 
ball in the same way as in Ref. [3], i.e. with vanishing 
initial transverse flow and with an initial energy density 
profile proportional to the transverse density of wounded 
nucleons, calculated from a Saxon- Woods nuclear den- 
sity profile with radius and surface thickness parameters 
Ro ^ 4.2 fm, C 0.596 fm for Cu and i?o = 6.37 fm, 
^ " 0.56 fm for Au nuclei. The energy density profile is 
normalized by a parameter eg = e(ro; r=&=0) giving the 
peak energy density in the center of the fireball for central 
collisions (impact parameter 6 = 0). cq is related to the 
peak wounded nucleon density in the same collisions by 
a factor k which is assumed to depend on energy but not 
on the size of the colliding nuclei. We here consider eg 
values that lead to final multiplicities covering the range 
accessible at RHIC and beyond, albeit perhaps not all 
the way to the Large Hadron Collider (LHC). 

As of now, the energy dependence of k cannot be cal- 
culated and must be determined empirically from the fi- 
nal charged hadron multiplicity dNch_/dy. Since dNch/dy 
counts the final entropy per unit of rapidity, including 
any entropy generated by viscous effects during the ex- 
pansion, the value of dNch/dy corresponding to a given 
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K will depend on the viscosity rj/s. We will see that 
the amount of entropy produced by viscous effects addi- 
tionally depends on system size, impact parameter and 
collision energy, but that all these dependences can, to 
good approximation, be absorbed in a single scaling func- 
tion, with parametric dependence on rj/s, that depends 
only on the multiplicity density (1/ S)dNch/dy: similar to 
V2l£, entropy production A<S/5o exhibits approximate 
"multiplicity scaling". However, this scaling function 
turns out to be non-linear. It therefore modifies the cen- 
trality dependence of the produced charged multiplicity, 
softening the observed increase with collision centrality 
of the produced charged multiplicity per pair of wounded 
Exploration of this important issue 



dN^i,/dy 



nucleons, 2 

requires an accurate modeling of the impact parameter 
dependence of the initial entropy density profile using, 
say, the Glauber or color glass condensate models. This 
is beyond the scope of the present article and will be left 
for a future study. 

Following the majority of previous studies [2-4, 45, 46], 
the viscous shear pressure tensor is initialized with its 
Navier-Stokes value tt™" = 2r]a^'^. Other initial con- 
ditions were studied in [1, 3], but the final observables 
were found to be insensitive to such variations [3]. The 
kinetic relaxation time for the kinetic evolution of the 
shear pressure tensor is taken as = c^r^"'*^ where 
^Boitz _ ^2 jg ^^f^ kinetic theory value for a classical 
gas of massless Boltzmann particles [38] and c^r is varied 
between j and 1. 

Decoupling from the hydrodynamic fluid is imple- 
mented by following the same procedure as described 
in [3]. We use the AZHYDRO algorithm [47] to find the 
freeze-out surface at constant decoupling temperature 
Idee = 130 MeV and calculate the final hadron spec- 
tra from the Cooper- Frye integral over this surface [48], 
with a distribution function that accounts for the remain- 
ing small deviations from local thermal equilibrium along 
that surface [3, 6]. Resonance decays are neglected, and 
only the elliptic flow of directly emitted pions is shown. 
To estimate the total charged hadron multiplicity, we 
take the directly emitted positive pions, multiply by 1.5 
to roughly account for multiplication by resonance decays 
at Tdec, then multiply by another factor 2 x 1.2 = 2.4 
to account for the negatives and roughly 20% of final 
charged hadrons that are not pions. A proper calculation 
of the resonance decay chain is computationally expen- 
sive and, for a systematic study like the one presented 
here that requires hundreds of runs of VISH2+1, beyond 
our presently available resources. 

Figure 1 shows the three equations of state (EOS) ex- 
plored in the present study. EOS I describes a non- 
interacting gas of massless particles, e = p/3. EOS Q is a 
frequently employed equation of state [25] that matches a 
non-interacting quark-gluon gas above Tc in a first-order 
transition (Maxwell construction) to a realistic hadron 
resonance gas (HRG) in chemical equilibrium below 
using a bag parameter B to adjust Tc to Tc = 164 MeV. 
SM-EOS Q is a shghtly smoothed version of EOS Q, see 




EOS I 
EOSQ 
SM-EOS Q 
EOSL 



5 6 
e (GeV/fm^) 

FIG. 1: (Color online) The equation of state. Panel (a) shows 
the pressure p as a function of energy density e and (in the 
inset) the squared speed of sound = || as a function of 

„2 



temperature T. Panel (b) shows 
density e. 



as a function of energy 



[3] for details. Since the discontinuities in the function 
c^(e) for EOS Q cause numerical problems in VISH2+1 
due to large velocity gradients near the interfaces be- 
tween QGP, mixed phase and HRG, we here use SM- 
EOS Q. 

EOS L matches the same hadron resonance gas be- 
low Tc smoothly in a rapid cross-over transition to lat- 
tice QCD data [49] above Tc- For the fit, the lattice 
data were plotted in the form p(e), interpolated and then 
smoothly joined to the p{e) curve of the HRG. As can 
be seen in the upper panel of Fig. 1 in the inset, our 
procedure is not fully thermodynamically consistent and 
leads to a somewhat different temperature dependence 
of cl below Tc than for EOS Q and SM-EOS Q. Since 
this only affects the flow dynamics below our decoupling 
temperature of Tdoc = 130 MeV, we have not put any 
effort into correcting this. For future comparisons of vis- 
cous hydrodynamic calculations with experimental data, 
the chemical equilibrium hadron resonance gas below Tc 
employed here must be replaced by a chemically non- 
equilibrated hadron gas whose particle ratios are frozen in 
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at the chemical decoupHng temperature Tchcm ~ T^, this 
has well-known consequences for the final hadron spec- 
tra and elliptic flow which can not be neglected [50] . We 
postpone this, together with a more careful and thcrmo- 
dynamically fully consistent matching to the lattice QCD 
data, to a future study. We note, however, that EOS L 
shown in Fig. 1 is quite similar to "EOS qp" studied in 
Ref. [51]. 



IV. EVOLUTION OF MOMENTUM 
ANISOTROPIES: SIMPLIFIED VS. FULL I-S 
EQUATIONS 

It has been previously observed that the results of 
Rcfs. [1] and [3] for the differential elliptic flow V2{pt), 
although both based on the Israel-Stewart 2"'^ order for- 
malism, seemed to disagree, our work [3] showing much 
stronger viscous suppression of V2 than that of P. & 
U. Romatschkc [1]. The resolution of this discrepancy 
was made difficult by the fact that the two groups not 
only used different versions of the Israel-Stewart equa- 
tion (3) as described in Sec. II, but also different initial 
conditions, different equations of state, and system sizes 
(Cu+Cu [3] vs. Au+Au [1]). In [3] we noted in a foot- 
note that the main reason for the observed differences 
seemed to be the different I-S equations used by the two 
groups. As we will see, this is only part of the story. In 
this sections we explore this question further and lead it 
to a complete resolution. 

Figure 2 shows the temporal evolution of the total mo- 
mentum amsotropy Sp = ]j.x:,_^rpyyi^ averaged over the 

transverse plane^ for two collision systems (Cu+Cu at 
b = 7fm on the left, Au-fAu at 6 = 7fm on the right) 
and three equations of state (EOS I (top), SM-EOS Q 
(middle), and EOS L (bottom)). The blue dashed lines 
at the top indicate the result from ideal fluid dynamics, 
the black and orange lines below show viscous fluid dy- 
namical results. The black lines show solutions of the 
full I-S equations, the orange ones for the simplified I- 
S approach; in each case several values of the kinetic 
relaxation time Ttt are explored. Note that our full I- 
S equations (3) do not use the identity (4) used in [1] 
which strictly holds only for conformal fluids (i.e. for the 
case of EOS I in Fig. 2). We have, however, tested the 
two expressions on the left and right of Eq. (4) against 
each other also for the other two equations of state (SM- 



^ Note that £p as defined here includes the effects from both fiow 
velocity and shear pressure anisotropies [3]. In Ref. [3] we de- 
noted it by e'p in order to distinguish it from the flow-induced 

momentum anisotropy . ° ^ Jj^i; which is based only on the 

(-'o +-'0 / 

ideal fluid part of the energy momentum tensor and neglects 
anisotropies in the local fluid rest frame caused by the shear 
pressure tensor tt™" . In the present work we drop the prime for 
convenience. 



EOS Q and EOS L) which are not conformally invariant, 
and found no discernible differences. Only for a very 
long relaxation time r^r = \2rj/ sT (not shown in Fig. 2) 
did we see for EOS L a difference larger than the line 
width, with our result for Sp lying slightly above the one 
obtained with the conformal approximation (4). 

Comparison of the black and red lines in Fig. 2 shows 
that the sensitivity of the momentum anisotropy Ep to 
the relaxation time Tjr is significantly larger for the sim- 
plified I-S equations (red) than for the full I-S equations, 
and that the T7r-dependence of Ep even has the oppo- 
site sign for the two sets of equations. With the full 
I-S equations, Ep moves slowly towards the ideal ffuid 
limit as r^r increases whereas with the simplified I-S equa- 
tions Ep moves away from the ideal fiuid limit, at a more 
rapid rate, resulting in a larger viscous suppression of the 
momentum anisotropy. In the limit r^r ^ 0, both for- 
mulations approach the same Navier-Stokes limit. The 
difference between full and simplified I-S equations is 
largest for EOS I which is the stiffest of the three stud- 
ied equations of state, causing the most rapid expansion 
of the fireball. For this EOS, the simplified I-S equa- 
tions allow for the largest excursions of tt™" away from 
its Navier-Stokes limit, causing a significant and strongly 
T7r-dcpendent increase of all viscous effects, including the 
suppression of the momentum anisotropy (Fig. 2) and el- 
liptic flow (see Fig. 4 below) and the amount of viscous 
entropy production (see Sec. VI). 

For the other two equations of state, SM-EOS Q and 
EOS L, the difference between full and simplified I-S dy- 
namics is much smaller, ranging from ~ 5% for Au+Au 
to ~ 15% for Cu+Cu for the largest value of Q-q/sT 
studied here. Note that the viscous suppression of Ep is 
much stronger for the smaller Cu+Cu collision system 
than for Au+Au. For SM-EOS Q and EOS L (which 
yield rather similar results for Ep, with differences not 
exceeding ~ 10%), the results from the full I-S equations 
(black lines) are almost completely independent of Ttt, 
even for the small Cu+Cu system. 

The insets in the two upper panels of Fig. 2 illustrate 
the different TTr-dependences for Ep in the full and sim- 
plified I-S formulations, by plotting the value of Ep for 
EOS I at a fixed time t — tq = 4dm/c as a function 
of T^. One sees that, for the investigated range of re- 
laxation times, the r^-dependence is linear, but that the 
slope has different signs for the full and simplified I-S 
equations and is much smaller for the full I-S system. 
Even though VISH2+1 cannot be run for much smaller 
values, due to numerical instabilities that develop as 
the Navier-Stokes limit = is approached, the lines 
corresponding to the full and the simplified I-S equations 
are seen to nicely extrapolate to the same Navier-Stokes 
point, as they should. For SM-EOS Q and EOS L, the 
corresponding lines may no longer be linear, due to phase 
transition effects, but are still characterized by opposite 
slopes for the simplified and full I-S approaches, with al- 
most vanishing slope in the full I-S case. This agrees with 
findings reported in [1, 52]. 
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Cu+Cu, b=7 fm, EOS I 




Au+Au, b=7 fm, EOS I 
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FIG. 2; (Color online) Time evolution of the total momentum anisotropy £p for two collision systems (left: Cu+Cu; right; 
Au+Au), three equations of state (top: EOS I; middle: SM-EOS Q; bottom: EOS L), and three values of the kinetic relaxation 
time r-n- as indicated (dotted, dashed and solid curves, respectively). The insets in the two top panels show the T,,r-dependence 
of the momentum anisotropy £p at fixed time r — ro = 4fm/c. See text for discussion. 



In Fig. 3 the effects of changing the system size, EOS, 
and form of 1-S equations on the differential elliptic flow 
V2{pt) for directly emitted pions is shown. The largest 
viscous suppression of elliptic flow (by almost 70% be- 
low the ideal fluid value at px — 2GeV/c) is seen for 



the small Cu+Cu system, evolved with SM-EOS Q and 
the simplified I-S equation. This is the result reported 
by us in [3]. The middle panel of Fig. 3 shows that this 
large V2 suppression is almost cut in half by going from 
Cu+Cu to Au+Au, the system studied in [1], even with- 
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FIG. 3: (Color online) Differential elliptic flow «2(pt) for directly emitted pions (i.e. without resonance decay contributions), 
comparing results for different collisions systems and equations of state, (a) Cu+Cu at 6 = 7fm with SM-EOS Q. (b) Au+Au 
at & = 7fm with SM-EOS Q. (c) Au-|-Au at b = 7fm with EOS L. Dashed (solid) lines correspond to ideal (viscous) fluid 
dynamics, with parameters as indicated. 



out modifying the EOS or the form of the I-S equations. 
Changing the EOS from SM-EOS Q [3] to EOS L (which 
is close to the one used in [1]) reduces the viscous V2 
suppression by another quarter, from about 40% to less 
than 30% below the ideal fluid limit at pt = 2. Finally, 
replacing the simplified I-S equations used in [3] by the 
full I-S equations employed by Romatschke [1] further 
reduces the suppression from about 28% below the ideal 
fluid to ^ 25% at p-p = 2GeV/c. This is consistent with 
the results obtained [1]. 

We conclude that the biggest contribution to the large 
difference between the results reported in Refs. [3] and [1] 
arises from the different collisions systems studied, with 
much larger viscous effects seen in the smaller Cu-|-Cu 
system than in Au-|-Au collisions. The next most im- 
portant sensitivity is to the EOS; for the most realistic 
EOS studied here, EOS L, the differences between using 
the full or simplified I-S equations with = Z-q/sT are 
only about 10% on a relative scale, or about 3% on the 
absolute scale set by the elliptic flow from ideal fluid dy- 
namics. For smaller t.^, this last difference would shrink 
even further. 

The sensitivity to details of the EOS documented by 
the middle and right panels of Fig. 3 gives an idea of 
how well one needs to know the EOS if one wants to ex- 
tract the speciflc shear viscosity rj/s from experimental 
data using viscous hydrodynamics. One might argue that 
the difference between a first order phase transition im- 
plemented through SM-EOS Q and a smooth crossover 
as in EOS L should be sufficiently extreme to cover the 
maximal theoretical uncertainty. In this case, Fig 3 tells 
us that the maximal theoretical uncertainty on the vis- 
cous suppression of V2 (and therefore on t]/ s) should be 
about 25-30%. This should be compared to the theoreti- 
cal error introduced by our present uncertainty of the ini- 
tial spatial source eccentricity e: e differs by about 30% 
between initializations based on the Glauber and color 
glass condensate (CGC) models [19, 52-54], resulting in 
a ~ 30% uncertainty of the total magnitude of the ellip- 
tic flow in ideal fluid dynamics. We further caution that 



recent discussions about the value of the critical tempera- 
ture for the quark-hadron transition [55, 56] introduce an 
additional moment of uncertainty which is perhaps not 
covered by the range between SM-EOS Q and EOS L 
studied here. Therefore, while we agree with the authors 
of Ref. [52] that the uncertainty about the initial source 
eccentricity dominates over uncertainties related to dif- 
ferent implementations of the I-S formalism, we think 
that the EOS should not be discounted prematurely as 
a possible source of significant additional theoretical un- 
certainty in the extraction of r//s. 



V. MULTIPLICITY SCALING OF v-i/e IN 
IDEAL AND VISCOUS HYDRODYNAMICS 

In this section we explore the multiplicity scaling (as 
defined in the Introduction) of the eccentricity-scaled el- 
liptic flow W2/e, comparing ideal fluid dynamics with that 
of near-minimally viscous fluids with specific shear vis- 
cosityf 



A. EOS I: conformal fluids with e = 3p 

We begin with the simple case of a conformal fiuid 
with the equation of state e = 3p (EOS I), without phase 
transition. In this case the speed of sound is a constant, 
independent of temperature T, = i. For the ideal 
fluid case, naive scaling arguments based on the scale in- 
variance of the ideal fluid equations of motion would thus 
predict a constant V2/e, independent of multiplicity den- 
sity {1/ S)dNch/dy. (The nuclear overlap area S is com- 
puted as = TTy^ (x^) (y^) where (...) denotes the energy 
density weighted average over the transverse plane.) The 
left panel in Fig. 4 clearly contradicts this expectation. 
Freeze-out at Tdoc = 130 MeV cuts the hydrodynamic 
evolution of the momentum anisotropy Ep short before 
the elliptic flow has fully saturated. As the left panel of 
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FIG. 4: (Color online) The eccentricity-scaled elliptic flow V2/e as a function of charged multiplicity density, {1/ S)dNch/dy, for 
a conformal fluid with EOS I. Results for Cu+Cu and Au+Au collisions with two different initial energy densities at a variety 
of impact parameters, as indicated in the legend, are superimposed. Results from ideal fluid dynamics (a) are compared with 
those from viscous hydrodynamics, using the full (b) and the simplified (c) Israel-Stewart equations, respectively. In all cases 
approximate, but not perfect multiplicity scaling is observed (see text for discussion). In panels (b) and (c), the ideal fluid 
results from the left panel are reproduced as gray symbols for comparison. 



the figure shows, this not only causes a strong suppres- 
sion of 1)2/^ at low multiplicity densities, where the time 
between beginning of the hydrodynamic expansion and 
freeze-out becomes very short, but it also breaks the mul- 
tiplicity scaling at high multiplicity density, albeit more 
weakly. At a fixed value of {l/S)dNch/dy, one sees larger 
V2/S for more central collisions initiated at lower collision 
energies (corresponding to smaller eo parameters) than 
for more peripheral collisions between the same nuclei at 
higher beam energies, and also for more central Cu-f Cu 
collisions (with a rounder shape) than for more peripheral 
Au-|-Au collisions (with a more deformed initial shape). 
We find that the larger V2/S values can be traced di- 
rectly to somewhat longer lifetimes of the corresponding 
fireballs, i.e. to the availability of more time to approach 
the saturation values of the momentum anisotropy and 
elliptic flow before reaching freeze-out. These freeze-out 
induced scaling violations in ideal fluid dynamics disap- 
pear at sufficiently high collision energies (i.e. large eg) 
where the momentum anisotropy has time to fully satu- 
rate in all collision systems and at all impact parameters, 
before freezing out. 

The middle and right panels in Fig. 4 show the ana- 
logous results for a minimally viscous fluid with 7 = 4^^ 
and kinetic relaxation time = ||;. Consistent with 
the discussion in the preceding section, the viscous sup- 
pression of the elliptic flow is seen to be stronger if the 
simplified I-S equations are used (right panel) than for 
the full I-S equation. (Although not shown, the curves 
in the right panel also show stronger sensitivity to the 
value of than those in the middle panel.) Along with 
the suppression of V2/S by shear viscosity we sec the ap- 
pearance of scale- breaking effects that increase in propor- 
tion to the overall suppression of elliptic flow: they are 
larger in the right than in the middle panel. Shear vis- 
cosity breaks the multiplicity scaling of V2 je because (as 
shown in the preceding section) viscous effects are larger 
in smaller collision fireballs. Consequently, if we compare 



different collision systems that produce the same charged 
multiplicity density (1/ S)dNch/ dy, we find smaller V2/S 
for Cu-f-Cu than for Au-I-Au collisions, and for periph- 
eral AuH-Au collisions at higher collision energy than for 
more central Au-f-Au collisions at lower collision energy. 

Viscous effects also generate entropy, i.e. they increase 
the final charged multiplicity dNch/dy. Comparing in the 
middle and right panels of Fig. 4 the gray (shaded) sym- 
bols from ideal fluid dynamics with the colored (solid) 
symbols for viscous hydrodynamics, points correspond- 
ing to the same collision system and impact parameter 
are seen to be shifted to the right. This enhances the 
scaling violations: for a given collision system, impact 
parameter and collision energy, viscosity decreases the 
eccentricity scaled elliptic flow V2/S, pushing the corre- 
sponding point downward in the diagram, and simulta- 
neously increases the entropy, pushing the corresponding 
point horizontally to the right. The combination of these 
two effects separates the curves for different collision sys- 
tems and energies farther in viscous hydrodynamics than 
in ideal fluid dynamics. 

B. Phase transition effects: EOS Q and EOS L 

Figure 5 shows the analogous results if the fluid evolves 
under the influence of an equation of state with a quark- 
hadron phase transition, EOS Q (top row) or EOS L 
(bottom row). Again approximate multiplicity scaling 
of V2/S is observed, but small scale-breaking effects are 
visible in both ideal and viscous hydrodynamics. For 
the equations of state with a phase transition, the scale- 
breaking effects are actually larger in the ideal than in the 
viscous case, i.e. in viscous hydrodynamics V2/£ shows 
better multiplicity scaling than in ideal fluid dynamics! 
We interpret the large scale-breaking effects in the ideal 
fluid case as a complication arising from interference be- 
tween the freeze-out process and the weak acceleration of 
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FIG. 5; (Color online) Same as Fig. 4, but for SM-EOS Q (top row) and EOS L (bottom row). For the ideal fluid case (a,d) 
an extended range of eo values up to eo = 120 GeV/fm'^ was studied, in order to show that V2/e eventually increases again at 
higher collision energies [25]. 



matter in the phase transition region. This interpretation 
is supported by a comparison between SM-EOS Q with 
its first-order phase transition (upper left panel in Fig. 5) 
and the smooth crossover transition in EOS L (lower left 
panel): for ideal fluids, the scale-breaking effects are ob- 
viously larger for SM-EOS Q than for EOS L. As already 
observed in [3] , shear viscosity effectively smears out the 
phase transition and reduces its effect on the dynamics. 
In Fig. 5 this is clearly seen on the left side of each panel 
(i.e. at small values of ^ ) where for the ideal fluid 
172 /e shows a non-monotonic peak structure [25] that is 
completely gone in the viscous case. 

Comparing Figs. 4 and 5, we see much smaller differ- 
ences between the full (middle panels) and simplified I-S 
equations (right panels) for SM-EOS Q and EOS L than 
for EOS I. This is consistent with our observations in 
Sec. IV where the largest differences between full and 
simplified TS equations was also seen for the rapidly 
evolving fireballs whose expansion is pushed by the very 
stiflt EOS I. 

It is interesting to observe that, for ideal fluids, EOS L 
leads to about 10% more elliptic flow under RHIC con- 
ditions than SM-EOS Q. The reason is that in the phase 
transition region EOS L is stiffer than SM-EOS Q. This 
plays an important role at RHIC because the softness of 
the EOS near Tc inhibits the buildup of elliptic flow ex- 
actly under RHIC conditions [25] . As a corollary we note 
that, if RHIC elliptic flow data exhaust ideal fluid pre- 



dictions made with SM-EOS Q [8], they will not exhaust 
ideal fluid predictions based on EOS L, thus leaving some 
room for shear viscous effects. 



C. Viscous suppression of D2: systematics 

Even at the highest collision energies (or eg values) 
studied in Figs. 4 and 5, the slope of V2/£ func- 
tion of ^^^^ remains positive, i.e. V2/£ continues to 
increase and evolve in direction of the asymptotic ideal 
fluid limit. This implies that at higher collision energies 
the importance of viscous effects decreases. This obser- 
vations parallels the one made in [3], namely that with 
increasing collision energy the pt range increases over 
which viscous hydrodynamic predictions for the single- 
particle momentum spectra can be trusted. The reason 
is in both cases that with increasing collision energy the 
time until freeze-out grows, and that (at least for con- 
stant rj/s as assumed here and in [3]) during the later 
stages of the expansion shear viscous effects are small. 

Figure 6 shows this more quantitatively. We plot the 
fractional decrease of the elliptic flow relative to its ideal 
fluid dynamical value, (4^^^=^^ - -i;™'=°"")/4'^<='^', as a func- 
tion of multiplicity density. Larger multiplicity densities 
lead to smaller viscous suppression effects. Larger vis- 
cosity results in stronger suppression of the elliptic flow. 
The suppression effects are weaker if the full I-S equa- 
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FIG. 6: (Color online) Viscous suppression of elliptic flow, 



function of {l/S)dNch/dy for 



Au+Au with EOS L, = Stj/sT and two values of 77/s as 
indicated. Solid (dashed) lines correspond to using the full 
(simplified) I-S equations, respectively. 



tions arc used than in the simplified approach of Ref. [3] 
(which, as discussed in Sec. IV, also suffers from strong 
sensitivity to r^). For minimal viscosity, rj/s = l/47r, 
the pT-integrated elliptic flow V2 in Au+Au collisions at 
RHIC is suppressed by about 20%. The suppression is 
larger at lower energies but will be less at the LHC. 



to confirm or disprove the theoretical predictions, but 
upcoming high-statistic runs at RHIC should be able to 
deliver them. 

Closer inspection of the two panels in Fig. 7 shows, 
however, that the theoretical scaling curves have the 
wrong slope: on the left side of the plot, i.e. for small 
multiplicity densities, the data seem to point towards 
larger specific shear viscosities 



> 3 X -r!- whereas 

47r 



on the right side of the plot, for -q 



1 dJVch 
S dy 



> 20 fm- 



thc experimental data require smaller shear viscosities, 
J < (1^2) X j^. But this is not at all unexpected: 
Collisions represented by points in the right half of the 
plot correspond to high collision energies and large ini- 
tial energy densities whose expanding fireballs spend the 
largest fraction of their life in the QGP phase. Fire- 
balls created in collisions represented by points in the 
left part of the diagram have smaller initial energy densi- 
ties and thus spend most of their time in the much more 
viscous liadronic phase [19]. A meaningful comparison 
between theory and experiment thus must necessarily ac- 
count for the temperature dependence of rj/s and its dra- 
matic increase during the quark-hadron phase transition 
[50]. This would lead to scaling curves in Fig. 7(b) with 
a larger slope that can better reproduce the data. What 
one can say already now is that the high-energy end of 
Fig. 7 requires very small specific shear viscosity rj/s for 
the QGP, of the same order as the minimal value postu- 
lated in [16, 17] (unless the initial source eccentricity e 
was strongly underestimated in the experimental data). 



D. A look at experimental data 



VI. MULTIPLICITY SCALING OF ENTROPY 
PRODUCTION IN VISCOUS HYDRODYNAMICS 



Figure 7(a) shows the famous experimental plot by 
Voloshin [21] with provides empirical evidence for mul- 
tiphcity scahng of V2/e. The lines labelled "HYDRO" 
are sketches for expectations from ideal fiuid dynam- 
ics, based on the calculations presented in [25] for V2 
in Au-|~Au collisions at fixed impact parameter 6 = 7fm 
as a function of multiplicity (parametrized by eo). They 
should be replaced by the curves shown in the left panels 
of Fig. 5. 

In Fig. 7(b) we present multiplicity scaling curves for 
V2/S obtained from viscous hydrodynamics with the full 
I-S equations. On a superficial level, the theoretical 
curves show qualitative similarity with the experimen- 
tal data, giving correct ball-park numbers if one assumes 
77/s ~ 0.24 ~ 3/47r. Interestingly, ignoring experimental 
error bars, one can see evidence for small scaling viola- 
tions in the experimental data whose pattern agrees with 
the theoretical predictions from viscous hydrodynamics 
(see discussion at the end of Sec. VA): the 62.5 AGeV 
Au4-Au data lie shghtly above the 200AGeV Au+Au 
points, and the 200 A GcV Cu+Cu points fall slightly be- 
low the 62.5 A GeV Au+Au data. Of course, these fine 
features of the experimental data are presently not sta- 
tistically significant; much more precise data are needed 



In the absence of shock waves, ideal fiuid dynamics 
conserves entropy. Correspondingly, the final multiplic- 
ity per unit rapidity is directly determined by the to- 
tal initial entropy per unit rapidity: 



dN 



const. X 



To / dx dy s{x,y;TQ). Numerical discretization of the 
hydrodynamic evolution equations introduces a small 
amount of "numerical viscosity" , however, which can not 
be fully avoided. To minimize numerical viscosity effects, 
the flux-corrected transport algorithm SHASTA used in 
the numerical solution of both the ideal and viscous fluid 
equations [3, 47] employs an "antidiffusion step" involv- 
ing a parameter called "antidiffusion constant" [57]. Nu- 
merical viscosity effects are maximal if this parameter is 
set to zero. In all our simulations we used 0.125 for the 
antidiffusion constant [57], resulting in about 0.3% en- 
tropy production by numerical viscosity in the ideal fluid 
case. In comparison with the 0(10—15%) entropy pro- 
duction in a viscous fluid with minimal shear viscosity 
(see below) this can be neglected'^. 



^ Effects from numerical viscosity depend on the spacing of the 
space-time grid used in the simulation. For the results presented 
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FIG. 7: (Color online) (a) The experimental observation of multiplicity scaling for V2/e, with data from Au+Au, Pb+Pb, and 
Cu+Cu collisions at RHIC, SPS and AGS [21]. (b) Theoretical prediction of approximate multiplicity scaling from viscous 
hydrodynamics using the full I-S equations, for three different (constant) specific entropy values rj/s = 0.08, 0.16, 0.24. 



Similar to the second paper of Ref. [46], we compute 
entropy production by exploiting the proportionality of 
final entropy to final charged multiplicity. We compute 



in this paper wo used Ax = Ay = 0.1 fm and At = 0.04 fm/c. 
To chock tho effects of numerical viscosity in the ideal fluid case 
wo also performed simulations with AZHYDRO where we either 
set the antidifi^usion constant to zero or increased and decreased 
the grid spacings [Ax, Ay, At) by a factor 2 to 4. To maximize 
numerical viscosity effects, we used EOS Q with a strong first or- 
der phase transition which generates shocks and associated large 
velocity gradients. We found that decreasing the grid spacing 
by a factor 2 has no visible effect on the average radial flow of 
the fluid in central Au+Au collisions but increases the momen- 
tum anisotropy of the ideal fluid in peripheral Au+Au collisions 
at fe = 7fm by 2—3%. A further reduction by another factor 2 
doesn't even affect the momentum anisotropy any more, indicat- 
ing that numerical viscosity effects have been basically reduced 
to zero. If we increase the grid spacing by a factor two, the ef- 
fects are a bit larger: the total entropy production by numerical 
viscosity increases by about 1%, the average radial flow at 6 = 
changes by 0.5%, and the momentum anisotropy experiences a 
relative suppression of about 5%. The largest effects are seen 
when maximizing the numerical viscosity at fixed grid spacing 
by setting the antidiffusion constant to zero. In this case the av- 
erage radial flow in central Au+Au collisions again changes only 
by 0.5% (which is negligible compared to the strong increase in 
transverse acceleration that we see in the viscous fluid), but the 
momentum anisotropy is reduced by about 10%, In the flnal pion 
V2 this is reflected by a py-depondent reduction that increases 
with pt, just like effects from real shear viscosity [1-4, 6]; the 
PT-integrated pion elliptic flow is reduced by less than 4% even 
in this extreme case. We conclude that numerical viscosity does 
not increase transverse acceleration but suppresses momentum 
anisotropy similar to real shear viscosity. For the parameters 
used in this paper, and for the equations of state studied here 
which do not have a sharp phase transition, numerical viscosity 
effects on the elliptic flow do not exceed 1—2% and can thus be 
neglected relative to effects from real shear viscosity. 



the final multiplicity dNch/dy for both ideal and viscous 
hydrodynamics and then equate the fractional increase 
in dNch/dy with the fractional increase in dS/dy. This 
ignores a small negative correction due to the viscous de- 
viation of the distribution function on the freeze-out hy- 
persurface from local equilibrium [38, 41] which slightly 
reduces the entropy per finally observed particle in the 
viscous case. The real entropy production is thus slightly 
smaller than calculated with our prescription. However, 
since on the freeze-out surface the viscous pressure com- 
ponents are small [3] , this correction should be negligible. 

We checked the above procedure by also directly in- 
tegrating the viscous entropy production rate 9 • s = 
TT^'^TT^i/ /277 over the space-time volume enclosed between 
the initial condition Cauchy surface and the final freeze- 
out surface. This method results in slightly larger en- 
tropy production, the relative difference amounting to 
about 0.7% (or about 0.07% in the absolute value of 
AiS/iSo) for central Au+Au collisions. Since the esti- 
mate from the final multiplicity gives a lower entropy 
production value even without accounting for the some- 
what smaller entropy per particle in the viscous case, we 
conclude that entropy production due to numerical vis- 
cosity must be a bit smaller in the viscous fluid than in 
the ideal one. This is not unreasonable, given the ob- 
servation in [3] that, compared to the ideal fluid case, 
the physical viscosity smoothens the strong velocity gra- 
dients near the quark-hadron phase transition, thereby 
presumably also reducing the effects of numerical viscos- 
ity. 

We note that our viscous evolution starts earlier (at 
To = 0.6 fm/c) than that of Ref. [46] (who use tq — 
Ifm/c). This earlier start results in larger entropy pro- 
duction fractions. As the inset in Fig. 8(b) shows, most of 
the entropy is produced during the early stage of the ex- 
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FIG. 8: (Color online) Entropy production AS, normalized by the initial entropy So, as a function of charged multiplicity 
density ^ . Calculations with VISH2+1 were performed for Au+Au and Cu+Cu collisions at various impact parameters 
and collision energies, using rj/s = 0.08, r,r = Srj/sT, and three different equations of state (EOS I, SM-EOS Q, and EOS L). 
(a) Simplified I-S equations, (b) Full I-S equations. The inset in panel (b) shows the entropy production as a function of time, 
for central Au+Au collisions with parameters as indicated in the legend. 



pansion. We have confirmed that the difference between 
Rcf [46] and the work; here is quantitatively reproduced 
by the entropy generated during the time interval from 
0.6 to i.Ofm/c, which can be calculated to excellent ap- 
proximation analytically [18] (using Eq. (D3) in Ref. [3]) 
by assuming boost-invariant longitudinal expansion with- 
out transverse flow during this period. 

Figure 8 shows the viscous entropy production AS, 
as a fraction of the initial entropy Sq, for Cu+Cu and 
Au+Au collisions at various impact parameters and col- 
lision energies, as a function of multiplicity density. One 
observes approximate multiplicity scaling of the frac- 
tional entropy production, with scaling functions that 
depend on the equation of state and, for non-zero ki- 
netic relaxation time, on the form of the Israel-Stewart 
equations used in the simulation. As with V2/S we see 
small scale-breaking effects, but generally the produced 
entropy fraction shows better multiplicity scaling than 
elliptic flow. The scale breaking effects for the viscous 
entropy production rate go in the same direction as with 
elliptic flow insofar as, at the same value of ^ , larger 
collision systems and more central collisions produce frac- 
tionally more entropy than smaller or more peripheral 
collisions, due to their longer lifetimes before freeze-out. 

For = 31] /sT, Figure 8 shows that the simplified 
TS equations (left panel) cause almost twice as much en- 
tropy production as the full TS system! Figure 9 clarifies 
that, when the simplified I-S equations are used, entropy 
production depends very sensitively on the kinetic relax- 
ation time Ttt, approaching the much smaller and almost 
completely r,r-independent entropy production rates of 
the full TS framework in the limit ^ 0. The large 
amount of extra entropy production for non-zero Ttt in 
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FIG. 9: (Color online) Sensitivity of the entropy production 
ratio AS /So shown in Fig. 8 to the kinetic relaxation time r-^, 
for the Au+Au collision system with eo — 30 GeV/fm (corre- 
sponding to a collision energy of 200 A GeV) . The three 
red curves (upper set) are for the simplified Israel-Stewart 
equations, the three black curves (lower set) for the full I-S 
equations. Solutions with the full I-S equations produce less 
entropy and show very little sensitivity to r-n. 



the simplified TS approach must thus be considered as 
unphysical. This is important because this artificial extra 
entropy production (caused by unphysically large excur- 
sions of the viscous shear pressure tensor tt™" away from 
its Navier-Stokes value tt"*" = 277(7™") manifests itself as 
additional charged hadron multiplicity in the observed 
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final state (seen as a shift of all points in the left panel 
of Fig. 8 towards larger values of -| ) . Since the final 
multiplicity is used to normalize the initial energy den- 
sity eo, this causes a significant distortion of the initial 
conditions corresponding to a given set of experimental 
data, affecting their physical interpretation. 



We conclude that using the full I-S equations is manda- 
tory if one wants to minimize artificial effects of shear 
viscosity on entropy production and elliptic flow in the 
realistic situation of non-zero kinetic relaxation times. 
(We note that, while the value of for the QGP created 
at RHIC is presently unknown, it can obviously not be 
zero.) From the right panel in Fig. 8 we see that this re- 
moves most of the large difference in entropy production 
between the rapidly exploding EOS I fireballs and their 
more leisurely expanding cousins that evolve under the 
influence of EOS Q or EOS L. StiU, even for the fuU I-S 
equations we see 15 — 25% differences between the entropy 
production rates for EOS Q (first order phase transition) 
and EOS L (rapid crossover transition). The differences 
are largest for the most central Au-|-Au and Cu-|-Cu colli- 
sions at top RHIC energies. The somewhat stiffer nature 
of EOS L near Tc causes the fireball to expand faster and 
with higher acceleration, leading to larger viscous effects 
than for EOS Q. The differences in entropy production 
caused by this variation of the EOS is of similar mag- 
nitude as its effect on the viscous suppression of V2/e 
discussed at the end of Sec. IV. 



An important comment relates to the negative overall 
slope of the scaling curves for entropy production shown 
in Figs. 8 and 9: Since peripheral collisions produce rela- 
tively more entropy than central collisions, and the pro- 
duced entropy is reflected in the final charged hadron 
multiplicity, the collision centrality dependence of hadron 
multiplicities is altered by viscous effects. When vis- 
cous effects are accounted for, the charged multiplicity 
dNch/dy will rise more slowly as a function of the number 
of participant nucleons A^part than for an ideal fluid with 
the same set of initial conditions. In a Glauber model 
parametrization of the initial conditions [8] this temper- 
ing effect will have to be compensated for by increasing 
the "hard" component in the initial entropy production, 
i.e. the component that scales with the density of bi- 
nary collisions and is thus responsible for the non-linear 
increase of dNch/dy with A'^part • In the color glass conden- 
sate approach [58] this non-linear rise is controlled by the 
centrality dependence of the saturation momentum scale 
Qs, with no free parameters to tune. It remains to be 
seen whether the success of the CGC model in describing 
the centrality dependence of dNch/dy [59] survives the 
inclusion of entropy (or multiplicity) producing effects 
resulting from shear viscosity during the evolution from 
the initial CGC to the finally observed state. 



VII. CONCLUDING REMARKS 

The main motivation for the work presented in this 
paper was provided by the experimentally observed mul- 
tiplicity scaling of the elliptic flow, shown in Fig. 7a, and 
its deviation at low multiplicities from ideal fluid dynam- 
ical predictions. We saw that many of the observed fea- 
tures arc qualitatively consistent with viscous hydrody- 
namic calculations as presented in this paper, and that 
the same calculations also predict approximate multiplic- 
ity scaling for viscous entropy production. Our studies 
revealed, however, that even for ideal fluid dynamics the 
multiplicity scaling of the elliptic flow is not perfect, with 
small scaling violations introduced by the freeze-out pro- 
cess which cuts the evolution of elliptic flow short. Even 
at RHIC energies, where the elliptic flow almost saturates 
before freeze-out, kinetic decoupling truncates the mo- 
mentum anisotropy at values slightly below their asymp- 
totic saturation value, and the deviations depend on the 
size of the colliding nuclei and the deformation of the flre- 
ball created in the collision through the time available for 
building elliptic flow before freeze-out. 

Shear viscosity strongly suppresses the build-up of mo- 
mentum anisotropy and elliptic flow, especially for low 
multiplicity densities, i.e. at large impact parameters, 
low collision energies or for small sizes of the colliding 
nuclei. This changes the slope of the multiplicity scal- 
ing curve for V2/S but preserves, to good approximation, 
its general scaling with ^ ^^^f^ ■ Violations of multiplicity 
scaling for V2/£ are somewhat larger for the viscous ex- 
pansion than for the ideal fluid (especially with EOS I), 
but remain small enough to be consistent, within statisti- 
cal errors, with the experimental observation of approxi- 
mate scaling. The slope of the approximate scaling curve 
and the spread around this curve caused by scaling viola- 
tion increase with the value of the speciflc shear viscosity 
r]/s and can thus be used to constrain it. 

Specifically, the observed scaling violations have the 
following features: At fixed multiplicity density 
viscous hydrodynamics predicts slightly larger elliptic 
flow V2/£ for larger collision systems or more central col- 
lisions than for smaller nuclei colliding at similar energy 
or more peripheral collisions between similar-size nuclei 
colliding at higher energy. Larger V2/£ values are associ- 
ated with longer lifetimes of the corresponding fireballs 
before freeze-out and thus also with larger relative en- 
tropy production. This correlates the scaling violations 
for V2/S observed in Figs. 5 and 7 with those for the 
relative entropy production AS/Sq seen in Fig. 8. The 
pattern of the predicted scaling violations shows qualita- 
tive agreement with experiment, although higher quality 
data are required to render this agreement statistically 
robust and quantitative. 

For a fixed (i.e. temperature independent) ratio rj/s, 
the slope of the multiplicity scaling curve for 1'2/e does 
not agree with experiment - the curves predicted by vis- 
cous hydrodynamics are too flat. The slope can be in- 
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creased by allowing 77/s to increase at lower tempera- 
tures: For small multiplicity densities (very peripheral 
collisions or low collision energies), the data seem to re- 
quire f > 3 X whereas at large multiplicity densities 
they appear to constrain the specific shear viscosity to 
values of ^ ^ x While this is qualitatively 

consistent with the idea that in high-multiplicity events 
the dynamics is dominated by the QGP phase (whose 
viscosity would thus have to be small, of order l/47r) 
whereas low-multiplicity events are predominantly con- 
trolled by hadron gas dynamics (which is highly viscous 
[19]), much additional work is needed to turn this ob- 
servation into quantitative constraints for the function 

The present study also resolves questions that arose 
from several recent publications of viscous hydrodynamic 
calculations which seemed to yield different results. We 
explored the effects of using different implementations 
of Israel-Stewart second order theory for causal rela- 
tivistic viscous hydrodynamics, by comparing the "sim- 
plified Israel-Stewart equations" previously used by us 
[3] with the "full Israel-Stewart equations" implemented 
by P. & U. Romatschke [1]. For the "simplified" ap- 
proach we found a strong sensitivity of physical observ- 
ables on the presently unknown kinetic relaxation time 

for the viscous shear pressure tensor tt™", in contrast 
to a much weaker and basically negligible T7r-dependence 
in the "full" approach. For non-zero the "simplified 
I-S equations" allow for large excursions of tt™" away 
from its Navier-Stokes limit tt™" ~ 2770-™". These ex- 
cursions are artificial and disappear in the Navier-Stokes 
limit — > which can, however, not be stably simulated 
numerically. They cause large viscous suppression effects 
for the elliptic flow and large amounts of extra entropy 
production (i.e. extra final hadron multiplicity). From 
our study we conclude that the "simplified I-S approach" 
should be avoided, and that a reliable extraction of rj/s 
from experimental data mandates the use of the "full 
Israel-Stewart equations" [1, 52]. (It is, however, per- 
missible to use the conformal fiuid approximation [1, 34] 
for the "full I-S equations" even if the fluid's EOS is not 
conformally invariant since the differences were found to 



be negligible.) 

In comparing our previous work [3] with that of others 
we also identified other factors that significantly influ- 
ence the creation of elliptic flow and thus help to ac- 
count for the observed differences. For a realistic equa- 
tion of state that implements a quark-hadron transition 
(here SM-EOS Q and EOS L), it turns out that a much 
more important effect than using the correct version of 
Israel-Stewart theory is the size of the colliding nuclei. 
At RHIC energies and for a realistic EOS, the viscous 
suppression effects for V2/e in Cu+Cu collisions are al- 
most twice as large as for the larger Au-t-Au collision 
system. Non-negligible differences in the amount of vis- 
cous V2 suppression arise also from details in the EOS, 
with a smooth crossover as implemented in EOS L giv- 
ing 25-30% less suppression than a first-order transition 
as in SM-EOS Q. Compared to system size effects and 
EOS uncertainties, the differences between "simplified" 
and "full" I-S theory are relatively small, affecting the 
viscous V2 suppression at the 10% level relative to each 
other. (The quoted percentages are for a fluid with min- 
imal viscosity rj/s = l/An and may be larger for higher 
viscosity.) The largest uncertainty, in any case, is con- 
tributed by our present lack of knowledge of the initial 
source eccentricity which contributes a theoretical error 
band of up to 30% on an absolute scale for V2 [19, 52-54]. 

We finally comment that the multiplicity dependence 
of viscous entropy production predicted by viscous hy- 
drodynamics (see Fig. 8) will modify the centrality de- 
pendence of charged hadron production. This issue will 
be studied more quantitatively in a forthcoming paper. 
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